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Abstract 

The single spin asymmetries for a longitudinally polarized lepton beam or a longitudinally polarized 
nucleon target in semi-inclusive deep-inelastic scattering are twist-3 observables. We study these 
asymmetries in a simple diquark spectator model of the nucleon. Analogous to the case of transverse 
qq , target polarization, non- vanishing asymmetries are generated by gluon exchange between the struck 

quark and the target system. It is pointed out that the coupling of the virtual photon to the diquark 
is needed in order to preserve electromagnetic gauge invariance at the twist-3 level. The calculation 

■ indicates that previous analyses of these observables are incomplete. 

o 

43 ' 1 Introduction 

a: 

Qh| The measurements of Ajjl (longitudinal target polarization) and Am (longitudinal lepton beam po- 
larization) by the HERMES [1, 2, 3] and CLAS [4] collaborations constitute the first clear evidence of 
non-vanishing single spin asymmetries (SSA) in semi-inclusive deep-inelastic scattering (SIDIS) off the 
nucleon. From a theoretical point of view, SSA in hard processes are very interesting because of their 
relation to time-reversal odd (T-odd) correlation functions (parton distributions and fragmentation 

■ functions). 

Already for more than a decade the existence of T-odd fragmentation functions is considered to be 
established [5]. In the meantime, explicit model calculations including final state interactions in the 
fragmentation process have provided non- vanishing results for such functions (see, e.g., Refs. [6, 7]). 
On the other hand, it has been shown only recently that nonzero T-odd parton distributions are 
compatible with time-reversal invariance of the strong interaction [8, 9] (see also Refs. [10, 11, 12] 
for related work). In DIS, T-odd parton distributions arise due to the exchange of longitudinally 
polarized gluons between the struck quark and the target system. This rescattering effect is encoded 
in the gauge link appearing in the definition of parton distributions. An alternative picture according 
to which T-odd parton distributions can be generated without rescattering of the struck quark [13] 
seems to be ruled out [14]. 

In particular, due to the recent developments [8, 9], the T-odd and transverse momentum depen- 
dent (fej_-dependent) functions f^p (Sivers function, describing the distribution of unpolarized quarks 
in a transversely polarized target) [15] and h\ (distribution of transversely polarized quarks in an 
unpolarized target) [16] may well exist. From a practical point of view both distributions can be 
considered as twist-2 functions, since they appear in observables at leading order of a 1/Q-expansion, 
where Q denotes the large scale of the hard process. For instance, fy^ enters the leading twist SSA 
Ajjt in SIDIS [16]. 
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Despite of the recent progress in understanding the nature of T-odd effects, a complete formalism 
(including subleading T-odd parton distributions) of T-odd twist-3 observables is still missing even at 
tree-level. So far, such effects have only been treated on the fragmentation side [17, 18]. This point 
may also be quite important for the description of the twist-3 asymmetries Ajjl and Aljj in SIDIS. 
With the exception of Refs. [19, 20], all present analyses/calculations of these observables are based 
on [18], i.e., they include only T-odd fragmentation functions (see e.g. Refs. [21, 22, 23, 24, 25, 26]). In 
this scenario one finds schematically Aljj ^ eH± and Ajjl cc Iil H^~, with the twist-3 T-even parton 
distributions e and Hl, and the twist-2 T-odd Collins fragmentation function H^- [5]. 1 

Analogous to the treatment of Ajjt presented in Ref. [8], we compute Am and Ajjl in the frame- 
work of a simple diquark spectator model of the nucleon, in order to investigate whether T-odd parton 
distributions may be relevant in these cases. The rescattering of the struck quark, which serves as 
the potential source of T-odd effects, is modelled by the exchange of an Abelian gauge boson. Such 
a study has already been performed in Ref. [27]. However, in [27] only Aljj has been computed ex- 
plicitly. Moreover, the calculation of [27] is not gauge invariant. To preserve electromagnetic current 
conservation also the coupling of the virtual photon to the diquark has to be included. 

Both Aljj and Ajjl turn out to be nonzero indicating that, in a factorized picture, T-odd distri- 
butions have to be taken into account. A first step in this direction has been made in [19], where it 
has been demonstrated that the T-odd distribution appears in the description of Aljj through a 
term hj~ E, where E is a twist-3 fragmentation function. As will be discussed below, our calculation 
of Aljj, however, cannot be identified with such a term suggesting that the formula of [19] for the 
beam SSA is not yet complete. 



2 Tree diagrams 

In order to study SIDIS off a spin-^ particle (for definiteness we think of a proton) in the framework 
of a spectator model we consider the process (compare also Ref. [8]) 

7*(?)+P(P,A)->g(pi,A') + s(p2). (1) 

In full SIDIS, both the quark and the spectator in the final state fragment into hadrons, where we 
are interested in the situation that one of the hadrons from the quark fragmentation is detected 
at low transverse momentum. However, for the study of possible T-odd effects related with parton 
distributions it is not necessary to include the fragmentation process in the calculation. We use the 
model of [8] with a scalar diquark spectator s. In this model the proton has no electromagnetic 
charge, and a charge e\ is assigned to the quark. The interaction between the proton, the quark and 
the spectator is described by a scalar vertex with the coupling constant g. 

We treat the process (1) in the Breit frame of the virtual photon. The proton has a large plus- 
momentum Q/x, where x = xjjj + 0{1/Q 2 ). The quark carries the large minus-momentum p\ m q~ 
and a soft transverse momentum Aj_. These requirements specify the kinematics: 

q=(-Q,Q,6 ± ), p=(9-, ^,0 ± ), (2) 

Q"' Q ' A± ) ' P2 ~ K~~^' Q(l-x) ' 

The expressions for q and p are exact, while for p\ and p2 just the leading terms have been listed. In 
particular, sometimes the 1/Q 2 corrections of p~{ and are needed which can be readily obtained 
from 4-momentum conservation. To simplify the calculation we consider massless quarks. 

1 In our calculation for Aul the target is polarized along the direction of its momentum (and the direction of the 
virtual photon). In experiments for the longitudinal target asymmetry, however, the polarization is along the direction 
of the incoming lepton. Both situations differ by a kinematical twist-3 term which is given by Aut- 
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Figure 1: Tree- level diagrams of the process in (1). In order to preserve gauge invariance both diagrams 
have to be considered. 



The tree- level diagrams of the process (1) are shown in Fig. 1. Their currents, depending on the 
helicities of the proton and the quark, read 



J (a,0)( A > A/ ) = e l3 



Oi - q) 2 



J (&,0)( A ' A ') 



-eig 



u(pi,\')'fWi-i)u(p, A), 

(2^-^)«( Pl ,A')«(p,A) 



(P2 - q) 2 - m 2 s 



(3) 
(4) 



We have defined the current by means of the scattering amplitude according to T = e^J 11 , with e 
denoting the polarization vector of the virtual photon. It is easy to check that current conservation 
holds for the sum of the two diagrams, i.e., 



If* \ J (a,0) + ^(6,0) 



0. 



(5) 



As long as one is just interested in leading twist observables, it is sufficient to consider diagram 
(a,0). This is for instance the case in the calculation of the transverse SSA Ajjt m Ref. [8]. The 
specific kinematics in Eq. (2) is the reason for the suppression of diagram (b,0) relative to (a,0). 
The propagator of the diquark in (b,0) behaves like 1/Q 2 , while there is no large momentum flow 
through the quark propagator in (a,0). Nevertheless, as we discuss in the following, for subleading 
twist observables diagram (b,0) can no longer be neglected. 

We compute the various components of the currents using the lightfront helicity spinors of [28]. 

For Jf n s one obtains 

O,o) 



J(a,0) ( A > A ') 
^(a,0)( A > A ') 

J (a,0)( A ' A ') 

J (a,0)( A ' A ) 



1 -x 



-eig 



Q 



Vx|A_l| A\ + m 2 L 
1 -x Q 



Mx^A 1 - i\A 2 y x ,\> - A (A 1 + iAA 2 )\_ A , 



-2ag 
-2ei 5 



|Aj_| A^ + fh? 
1-x |A i I 



\fx A\ + fh? 
1-x |A i I 



Mx (iXA 1 + A 2 )5 XjX , -i(- iXA 1 + A 2 ) 5 x ,-\> 
Mx 8 XX i — f AA 1 + z A 2 ) 8 X X i 



y/x A\ + m 2 
with rh 2 = x(l — x) ( 



Mx5 



A, A' 



^(i-4 7 ±4Vaa1 + ,A 2 V 
A 2 i V M 2 (l-x)A / 



A,-A' 



(6) 
(7) 

(8) 

,(9) 



M + 



inz 



1 — x 



3 




Figure 2: One-loop diagrams for the imaginary part of the process in (1). In each diagram the possible 
on-shell intermediate state is indicated by a thin dotted line. The interaction between the quark and 
the spectator is modelled by the exchange of an Abelian gauge boson. 



One observes here the well-known result that for DIS off a spin-| particle the transverse current is 
dominating in the Breit frame. For the second tree-graph we find 

J+ (A,Y) = -e lg 2 ~ * (A A 1 +iA 2 }6\-\i , (10) 
' V x I A j. | v ' 



<W A ' A ') = eiSj^Pj^ + iA 3 )^.*. (11) 



The transverse components J) b0 \ are proportional to l/Q and, hence, indeed are suppressed by a 
factor l/Q 2 compared to J? ^ as expected. Therefore, these terms are not relevant for the discussion 
of twist-3 observables. In contrast, the plus- and the minus-components for both diagrams are of 
the same order. In this case, the suppression of (b,0) caused by the propagator of the diquark is 
compensated by a factor Q at the photon-diquark vertex and the fact that u(pi, A') u(p, A) oc Q. Even 
though diagram (b,0) is not compatible with the parton model, since large momentum transfers at 
the proton-quark-diquark vertex are allowed, a consistent calculation of twist-3 observables in the 
spectator model must contain this contribution. We note that our results obey the gauge invariance 
constraint 

J {a,0) + ^(6,0) = J {a,0) + ^(6,0) ■ ( 12 ) 
Including by hand a formfactor at the proton-quark-diquark vertex in order to suppress large momen- 
tum transfers destroys the gauge invariance. 



3 One-loop diagrams 

To obtain nonzero SSA in a spectator model one has to go beyond the tree-level approximation and take 
the rescattering of the quark into account. For our purpose it is sufficient to model this effect by one- 
photon exchange. Since the imaginary part of the one-loop amplitude is needed for the computation of 
SSA, just the two diagrams in Fig. 2 have to be considered. Self-energy and vertex correction diagrams 
are relevant for the real part of the amplitude, but cannot acquire an imaginary part, because we are 
dealing with either on-shell or even space-like (internal) lines. In order to avoid infrared singularities 
at intermediate steps of the calculation we assign a mass fx to the photon. The final results for Am 
and Ajjl must be infrared-finite which serves as a non-trivial check of the calculation. The currents 



4 



of the diagrams in Fig. 2 are given by 

d 4 k 



u(pi, A') 2jfc- JS 7" - rf) u(p, A) 



Jf M) (A,A') = -i( ei fgj 



[k 2 + ie] p - g) 2 + ie] [{p + q- k) 2 - m 2 s + ie] [(k - Pl ) 2 - n 2 + ie] ' 



(2tt)< 



(13) 



(14) 



u(pi, A') (^1+ 2^2- #) # (2p" + - 2*^) u(p, A) 



[fc 2 + ie] [(p — /c) 2 — m 2 + ie] [{p + q— k) 2 — m 2 + ie] [(k — pi) 2 — fi 2 + ie] 
Applying Cutkosky-rules to calculate the imaginary part one can verify the gauge invariance condition 

^( imj M)+ imj rM))=°- ( i5 ) 

It is obvious that the full current (including the real part) for the sum of both diagrams is not gauge 
invariant. 

The calculation of the imaginary parts has been performed similar to the study of T-odd fragmen- 
tation in Ref. [29]. We refrain from giving any details and just quote the final results. For diagram 
(a,l) we find 



1^4,1) (A, A') 



ImJ 2 ajl) (A,A') 



(eifg 1-x 



8vr ^/x\A ± 
1 



Q 



^Mx^-iAA 2 )^, 



(16) 



2A\ 



(Ai " - 2 + M 2 ) + In ^ j a(a* + *AA 2 ) \_ x , 



(e^g 1-x 



8vr 



(17) 



2A\ 



(A3 



fa + n' 



InL 

7x 



+ ki- 



rn 



ImJ+ al) (\,\') = - 



{eifg 1-x 



Mx [ ( A\ - m 2 + fj 2 



InL 



+ ln- 



iAA 1 + A 2 j 5 x ,_y 

~ 2 - 
m 



I 1 



+ l) (AA^A 2 )^, 

, , A=(Ai + m 2 -^ 2 ) 2 + 4Ai / , 2 

A 2 | + m 2 + u 2 - yA v ' 



(18) 



_ 2 lnL Q 2 (l-x) 

+2( mr—= - In 5 

VA /U z x 

A^ + m 2 + /i 2 + \/A 



with L = — 



The plus-component of diagram (b,l) is given by 



ImJ+ 1) (A,A') = - 



{eifg 1 



8vr v^|A ± | 



(2 - x) In ^ 

fl Z X 



2(l-x))(AA 1 + i A 2 )5 A ,_ 



A' • 



(19) 



For the one-loop calculation we make use of gauge invariance to eliminate the minus-component of 
the current. The Q-behaviour of the one-loop expressions corresponds to the one of the tree-graphs. 
Note also that the plus-component of the currents for both diagrams contains a In Q 2 -term, which is 
not compatible with the parton model. From our results for the transverse currents in Eqs. (16,17) 
we were able to reproduce the transverse target SSA computed in Ref. [8] (up to an overall sign). 
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4 Spin asymmetries 



Eventually, we proceed to the calculation of Am and Ajjl- The full cross section in DIS (including 
the leptons) in the one-photon exchange approximation can be expressed in the standard form 



with the lepton tensor 



a oc L. U W^ , 



IT = 2{lH w + n'" -^g^ + i\ e e^ pa q p l 



(20) 



(21) 



In Eq. (21), the 4-momentum of the incoming (outgoing) lepton is denoted by I (I'), and Q 2 = — (I — I') 2 . 
The hadron tensor is obtained from the above currents according to 



(>)V. 



(22) 



Now, we exploit gauge invariance of both the lepton and hadron tensor, take q from Eq. (2), choose the 
lepton momenta to be in the xz— plane, and ignore a contribution of iy ++ which for our calculation 
is at least suppressed by a factor l/Q 2 relative to the leading term in the cross section. This allows 
us to write the cross section as 



oc 



4Q 2 



2 

- (2 - y ) VT^yW+ 1 + {i- y + y -)wl 1 + 

+i\Jy^/T^Wf-y(l- 



v —wf 

4 b 



(23) 



where Wg" = {W^ + W v ^)/2 and W%" = (W^ - W vtl )/2 represent the symmetric and the antisym- 
metric part of the hadronic tensor respectively. We also used the standard definition y = p ■ q/p • I. 
The first line in Eq. (23) is relevant for the target spin asymmetry, and the second one for the beam 
asymmetry. Actually, it turns out that the purely transverse components of the hadronic tensor 
(Wg 1 , Wcj 2 , W^ 2 ) don't contribute to the spin-asymmetries at twist-3 level. 
In order to specify the asymmetries we define 



W, 



n 

uu 



wi 2 LU 



Z A, A' 

7 E f( J+ ( A ' A'))^ 2 (A, A') - (J 2 (A, X')) ] J + (A, A') 

A, A' L 

iEf((^(T.A / ))V(T,v)-(j + a,v)) t J l a > A' 
+((j i (T,A')) t j+(T,A')-(j i a,A')) t j+a,A' 



(24) 
(25) 
(26) 



where polarization "|" in (26) means polarization along the positive z-axis, i.e., along the direction of 
the target momentum. The element Wyjj is given by the tree- level result for diagram (a,0) in Eq. (6). 
Nonzero contributions to Wj[ 2 LU and Wg^ L are generated by interference of the tree-level amplitude 
with the imaginary part of the one-loop amplitude. While the transverse currents in (25) and (26) 
are obtained from diagrams (a,0) and (a,l), all four diagrams contribute to the plus-component of the 
current. The final results for the asymmetries read 2 



A 



LU 



(eQ 2 Aj + rh 2 A 2 
4vr M 2 x 2 + A\Q 



M 2 x 2 



m 



2(1 -x) 



In 



A 2 , +fn 2 



(27) 



2 As a reference we mention that the transverse target spin asymmetry of [8] is given by Aut = W^/W^jj, with 



W S UT defined analogous to Eq. (26). 



6 



A 



UL 



Wbu 



(ei) 2 K\ + m 2 A 2 
4vr M 2 x 2 + A 2 ~Q 



x 



In 



Q 2 (i 



2(1 -x) (A^ + m 2 ; 



1 

I 7 



mV 



to 2 



2-z 



2(1 



In 



2(1 

Q 2 (i- 



(A^ + m 2 )x 



In 



A 2 ; +m 2 



m 2 



(28) 



We would like to add some remarks: 



• An explicit nonzero result for Alu in the framework of the diquark spectator model has already 
been obtained in Ref. [27]. Our calculation shows that both Am and Aul remain finite once all 
diagrams required by electromagnetic gauge invariance are taken into consideration. 

• We believe that the effect which generates Alu m our calculation is not related to a term 
proportional to h^- E discussed in Ref. [19]. While hj~ is chirally odd, we have summed over the 
polarizations of the outgoing quark. 

• Since the asymmetries are proportional to A 2 = | Aj_ | sin ^ we expect in full SIDIS an effect 
proportional to sin 4>hi where <ph is the azimuthal angle of the produced hadron. The mechanisms 
which have been discussed so far in the literature in connection with Alu and Ajjl [17, 18, 19] 
show the same sin 0/j-behaviour. In addition, the different contributions to the asymmetries have 
the same y-dependence (y-^/1 — y for Alu an d (2 — y)y/l — y for Aul, see Eq. (23)). 

• In the final results for the asymmetries we have performed the limit fx — > without encountering 
a divergence. We agree with the observation made in Ref. [27], that the contribution from 
diagrams (a,0) and (a,l) to Alu is separately infrared-finite. This behaviour, which holds for 
Aul as well, seems to be accidental. 



5 Summary and conclusions 

In summary, we have calculated the twist-3 single spin asymmetries Alu and Aul for semi-inclusive 
DIS off a nucleon target in the framework of a simple diquark spectator model. Our study completes 
the previous work in Ref. [27], where only Alu has been computed explicitly. Moreover, the treatment 
in [27] is lacking electromagnetic gauge invariance. 

Both Alu and Aul turn out to be nonzero. Although the spectator model calculation contains 
contributions which are not compatible with the parton model, the non- vanishing results indicate that 
T-odd distributions have to be included in a factorized description of the asymmetries. So far this has 
only been done partly in the literature. In fact, we have argued that apparently none of the present 
analyses/calculations of Alu and Aul within the parton model is complete. Mainly for two reasons 
we feel confident to make such a speculation: first, within our calculation non-zero asymmetries arise 
already from the diagram (a,l), whose kinematics is compatible with the parton model. Second, 
there is no reason why the asymmetries should not contain higher order T-odd distribution functions. 
The status of the parton model formulae for Alu and Aul needs to be clarified before any definite 
conclusion can be extracted from the data. 



Note added: After this work has been completed, a revised parton model analysis for Alu and 
Aul appeared [30]. The analysis confirms our suspicion that both asymmetries should contain an 
additional term with a twist-3 T-odd distribution function, which have not been taken into account 
in the literature before. 



7 



Acknowledgements: We are grateful to J.C. Collins and N. Kivel for discussions. The work has been 
partly supported by the Sofia Kovalevskaya Programme of the Alexander von Humboldt Foundation, 
the DFG and the COSY-Juelich project. 

References 

[1] HERMES Collaboration, A. Airapetian, et al., Phys. Rev. Lett. 84, 4047 (2000). 

[2] HERMES Collaboration, A. Airapetian, et al., Phys. Rev. D 64, 097101 (2001). 

[3] HERMES Collaboration, A. Airapetian, et al., Phys. Lett. B 562, 182 (2003). 

[4] CLAS Collaboration, H. Avakian, et al., Phys. Rev. D 69, 112004 (2004). 

[5] J. C. Collins, Nucl. Phys. B396, 161 (1993). 

[6] J. C. Collins, G. A. Ladinsky, hep-ph/9411444. 

[7] A. Bacchetta, R. Kundu, A. Metz, P. J. Mulders, Phys. Lett. B 506, 155 (2001). 

[8] S. J. Brodsky, D. S. Hwang, I. Schmidt, Phys. Lett. B 530, 99 (2002). 

[9] J. C. Collins, Phys. Lett. B 536, 43 (2002). 

[10] S. J. Brodsky, D. S. Hwang, I. Schmidt, Nucl. Phys. B642, 344 (2002). 

[11] A. V. Belitsky, X. Ji, F. Yuan, Nucl. Phys B656, 165 (2003). 

[12] D. Boer, P. J. Mulders, F. Pijlman, Nucl. Phys. B667, 201 (2003). 

[13] M. Anselmino, V. Barone, A. Drago, F. Murgia, hep-ph/0209073. 

[14] P. V. Pobylitsa, hep-ph/02 12027. 

[15] D. W. Sivers, Phys. Rev. D 41, 83 (1990); 
Phys. Rev. D 43, 261 (1991). 

[16] D. Boer, P. J. Mulders, Phys. Rev. D 57, 5780 (1998). 

[17] J. Levelt, P. J. Mulders, Phys. Lett. B 338, 357 (1994). 

[18] P. J. Mulders, R. D. Tangerman, Nucl. Phys. B461, 197 (1996); 
Nucl. Phys. B484, 538 (1997) (E). 

[19] F. Yuan, Phys. Lett. B 589, 28 (2004). 

[20] L. P. Gamberg, D. S. Hwang, K. A. Oganessyan, Phys. Lett. B 584, 276 (2004). 

[21] K. A. Oganessyan, H. R. Avakian, N. Bianchi, A. M. Kotzinian, hep-ph/9808368. 

[22] A. V. Efremov, K. Goeke, M. V. Polyakov, D. Urbano, Phys. Lett. B 478, 94 (2000). 

[23] E. De Sanctis, W. D. Nowak, K. A. Oganessyan, Phys. Lett. B 483, 69 (2000). 

[24] B. Q. Ma, I. Schmidt, J. J. Yang, Phys. Rev. D 63, 037501 (2001). 

[25] A. V. Efremov, K. Goeke, P. Schweitzer, Phys. Lett. B 522, 37 (2001); 
Phys. Lett. B 544, 389 (2002) (E). 



8 



[26] A. V. Efremov, K. Goeke, P. Schweitzer, Phys. Rev. D 67, 114014 (2003). 

[27] A. Afanasev, C. E. Carlson, hep-ph/0308163. 

[28] G. P. Lepage, S. J. Brodsky, Phys. Rev. D 22, 2157 (1980). 

[29] A. Metz, Phys. Lett. B 549, 139 (2002). 

[30] A. Bacchetta, P. J. Mulders, F. Pijlman, Phys. Lett. B 595, 309 (2004). 



9 



